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We discuss observational constrains coming from supernovae imposed on the behaviour of the 
Randall- Sundrum models. We test the models using the Perlmutter SN la data as well as the new 
Knop and Tonry/Barris samples. The data indicates that, under the assumption that we admit 
zero pressure dust matter on the brane, the cosmological constant is still needed to explain current 
observations. We estimate the model parameters using the best-fitting procedure and the likelihood 
method. The observations from supernovae give a large value of the density parameter for brane 
matter Q\,o — 0.01 as the best fit. For high redshifts z > 1.2, the difference between the brane 
model and the ACDM (Perlmutter) model becomes detectable observationally. From the maximum 
likelihood method we obtained the favored value of Qx.o ~ 0.004±0.016 for fi^o = and f2 m ,o = 0.3. 
This gives the limit il\ t o < 0.02 at la level. While the model with brane effects is preferred by 
the supernovae type la data, the model without brane fluid is still statistically admissible. We also 
discuss how fit depends on restrictions of the sample, especially with respect to redshift criteria. 
We also pointed out the property of sensitive dependence of results with respect to choice of M 
parameter. For comparison the limit on brane effects which comes from CMB anisotropics and 
BBN is also obtained. The uncertainty in the location of the first peak gives a stronger limit 
fi A ,o < 1-0 • 10~ 12 , whereas from BBN we obtain that Qa,o < 1-0 ■ 10 27 . However, both very 
strict limits are obtained with the assumption that brane effects do not change the physics in the 
pre-recombination era, while the SN la limit is model independent. 

We demonstrate that the fit to supernovae data can also be obtained if we admit the phantom 
matter p — — (4/3) £ on the brane, where this matter mimics the influence of the cosmological 
constant. We show that phantom matter enlarges the age of the universe on the brane which is 
demanded in cosmology. Finally, we propose to check for dark radiation and brane tension by the 
application of the angular diameter of galaxies minimum value test. 

PACS numbers: 04.20. Jb,04.65.+e,98.80.Hw 



I. INTRODUCTION 



The idea of brane universes has originated from Hof ava 
and Witten Q followed by Randall and Sundrum @] • In 
the Randall-Sundrum scenarios of the brane-world cos- 
mology the large extra dimensions can solve the mass 
hierarchy problem of the standard model. In this model 
the basic cosmological equations are modified by the pres- 
ence of some extra terms which are derived from the fact 
that our universe is treated as a three-brane - to which all 
the gauge interactions are confined - embedded in a five- 
dimensional, anti-de Sitter space, in (the whole of) which 
only gravity can propagate. The cosmological evolution 
of such brane universes has been extensively investigated 
01 01 IE @ (f° r a review see 0)- For example, the 
issues related to the cosmological constant problem 0, 
the cosmological perturbations [E El ! inflationary solu- 
tions homogeneity, and flatness problems [2> E3 
were discussed. 
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Brane models admit new parameters which are not 
present in standard cosmology (brane tension A and dark 
radiation U). From the astronomical observations of su- 
pernovae la Hl[ll E [H E 13 one knows that the 
universe is now accelerating and the best-fit model is for 
the 4-dimensional cosmological constant density param- 
eter f2A 4 ,o = 0.72 and for the dust density parameter 
^m,o = 0-28 (index "0" refers to the present moment of 
time). In other words, only the exotic (negative pressure) 
matter in standard cosmology can lead to this global ef- 
fect. On the other hand, in brane models the g 2 quadratic 
contribution in the energy density g even for a small neg- 
ative pressure, contributes effectively as the positive pres- 
sure, and makes brane models less accelerating. In this 
paper we argue that in order to avoid this problem one 
requires much stronger negative pressure p < — p (p han- 
tom) matter to be present on the brane (cf. Ref. |2Cil2"H ). 

We concentrate on some observational constraints on 
the g 2 term, which depending on the type of matter 
present scales with the cosmic scale factor as a~ 67 , 
where 7 is the barotropic index in the equation of state 
p = (7 — l)g (p - the pressure and g - the energy den- 
sity). In fact, the g 2 term scales as a~ 6 (similarly as 
stiff-fluid in standard cosmology) for dust 7=1 mat- 
ter, as a~ 8 for radiation 7 = 4/3 and as a 2 for phantom 
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matter 7 = —1/3. 

In this paper we discuss the observational constraints 
imposed on the brane-world cosmologies by supernovae 
type la observations. Preliminary results of such analysis 
were made in |32j . Since we do not know how the exotic 
physics of brane theory works in the early stages of the 
Universe (for example, during structure formation), then 
all estimations of the brane influence onto the Cosmic Mi- 
crowave Background (CMB) and Big Bang Nucleosynthe- 
sis (BBN) are based on the assumption that this exotic 
physics does not change physical processes afterwards. 
The advantage of SNIa data is its independence of the 
physical processes in the early universe. Therefore, de- 
spite the fact that weaker limits from SNIa observations 
can be obtained, they are in this sense more valuable. 
We demonstrate that the brane model fits the SN la data 
very well. We also show that these observations indicate 
a substantial contribution from the brane, although the 
observational constraints from BBN and CMB restrict 
this contribution to a much higher degree. However, to 
derive the limits from BBN and CMB we assume that 
brane effects did not change physics before the recom- 
bination epoch. On the other hand, the less restrictive 
limit we obtained from SN la data is not so strongly de- 
pendent on the assumptions of a model. 

We argue that the admission of the brane effects for 
the dust matter on the brane does not fit supernovae 
data without the presence of the cosmological constant, 
but we can avoid this problem if we admit phantom mat- 
ter on the brane. Also, we argue that in the near future, 
when new supernovae data is available (SNAP3), the hy- 
pothesis of important brane effects at present (f2A,o > 0) 
will be tested. 

In Section II we present the basic equations of brane 
cosmology which are studied dynamically in Section III. 
Sections IV and V are devoted to a redshift-magnitude 
relation for brane universes and its comparison with su- 
pernovae type la data. In Section VI we discuss the angu- 
lar diameter size of a galaxy for brane universes together 
with the age of the universe problem. In Section VII we 
obtain the restrictions onto the brane universes from the 
observations of the Doppler peaks in CMB while in Sec- 
tion VIII we give the restrictions which come from BBN. 
In Section IX we give our conclusions. 



models is obtained from the 5-dimensional Einstein equa- 
tions with brane location at y — as 

^G^ = 4l-4 5) 9^ + S(y)(-X {i) V + (4) V)] , (1) 

where ^G^ is the 5-dimensional Einstein tensor, ^g^v 
is the 5-dimensional metric, y is an orbifold coordi- 
nate, A5 is the 5-dimensional cosmological constant, 
= ^g^v — n^riy is 4-dimensional induced met- 
ric, n a is the unit vector normal to the brane, ^T^ v 
is the 4-dimensional energy-momentum tensor, A is the 
brane tension, and n\ = 8-k^Gn = Stt/^M^ is the 
5-dimcnsional Einstein constant. The main point is 
that the 5-dimensional Planck mass can be much less 
from the 4-dimensional Planck mass ^M p <C ^ M p = 
1.2 • 10 19 GeV and so the electroweak scale (5) M P ~ TeV 
may be reached in the Earth accelerators. 

The 4-dimensional Einstein equations on the brane are 
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where A4 the 4-dimensional cosmological constant and 
k 2 {4) = 8nG = 8tt x lCr 38 (GeV)- 2 is the 4-dimensional 
Einstein constant. For the homogeneous and isotropic 
Friedmann models the Einstein field equations reduce to 
a generalized Raychaudhuri equation [j| 
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-H 2 



p + 3p + (2p + 3p) 



A4 
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where the dark radiation U results from the Weyl tensor 
contribution to 5-dimcnsional geometry which is related 
to the fact of non- vanishing of the energy-momentum ten- 
sor of matter on the brane. The dark radiation obeys the 
conservation law 



II. BASIC EQUATIONS OF BRANE 
COSMOLOGY 



U = -AHU 



(7) 



The conservation law for matter in the form of perfect 
fluid is 



Hofava and Witten [lj and Randall and Sundrum |2j 
provided us with an alternative mechanism of compact- 
ification of extra dimensions which is different from the 
conventional Kaluza-Klein scheme because extra dimen- 
sions are segments (called orbifolds) rather than compact 
spaces like the circle. Based on the Randall-Sundrum 
brane model we can derive the effective 4-dimensional 
Einstein equations by projecting 5-dimensional metric 
onto the brane world-volume. The dynamics of brane 



p = -3(p + p)H 
which for barotropic fluid gives 

p = -3-fHp , 



(8) 



(9) 



and shows that the conservation law for equation JJJ) is 
simply the conservation law for radiation (7 — 4/3), ex- 
cept U can take on both positive and negative values. 
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In Ref. [23j we gave the formalism to express dynami- 
cal equations in terms of dimensionless observational den- 
sity parameters fi. Following Refs. j^ESES we intro- 
duce the notation useful for this purpose. In this nota- 
tion the Friedmann equation for brane universes takes 
the form which is the first integral of Equations JHJ) 



K 2 K 2 
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(10) 



(4) 



where a(t) is the scale factor, k = 0,±1 the curvature 
index. After imposing conservation law we have 



ga 3j = const. 



(11) 



which brings the Eq. i)l(J)) to the form 



1 fda\ _ C G r C x k A 4 C u (m 
a 2 \dt) a 3 f + a 67 a 2 3 a 4 ' 1 ' 

In Eq. 1)121) we have defined the appropriate constants 
C GR = (l/3) K 2 4) a^g, C x = 1/6A • ^ 4) a 6 V, C u = 

2/ K( 4 )A- a 4 l4 , and Cgr is a of general relativistic nature, 
C\ comes as contribution from brane tension A, and Cu 
as a contribution from dark radiation. Though in Refs. 
[23l |2tJ the Eq. (|T2|l was studied using qualitative meth- 
ods we briefly discuss here the cases 7 = (cosmological 
constant), 7 = 1/3 (domain walls) and 7 = 2/3 (cos- 
mic strings) which can exactly be integrable in terms of 
elementary or elliptic [24[ functions. 

The first case 7 = is the easiest, since the first two 
terms on the right-hand-side of (|12fl play the role of cos- 
mological constants similar to A4 



A_ 4 
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The next two cases involve terms which were already in- 
tegrated in the context of general relativity. For 7=1/3 
(domain walls on the brane) the general relativistic term 
Cgr. in 1)121) scales as domain walls in general relativity 
while the term with C\ scales as cosmic strings (curva- 
ture) in general relativity, i.e., 



1 (da 
a 2 \~dt 



C\ — k Cgr. A4 Cu .v 

2 o 4 ' V-^J 



For 7 = 2/3 (cosmic strings) the general relativistic term 
Cgr m H12I) scales as cosmic strings in general relativity, 
while the term with C\ scales as radiation in general 
relativity (compare |28|,|44(), i- e - ; 



1 ( da 



Cu + C\ C G r - k A4 
a 4 a 2 3 ' 



(15) 



with an effective curvature index k' = k — Cgr- Then, 
the problem of writing down exact solutions, which are 
elementary, reduces to the repetition of the discussion of 



Ref. |2J|. We will not be doing this here. For other 
values of 7 = 4/3; 1; 2 the terms of the type 1/a 8 and 
I /a 12 appear, and the integration involves hyperelliptic 
integrals. 

Let us also consider the case of 7 = —1/3 (phantom) 
H3, i.e., 

1 ( da\ 2 2 k C u A 4 

=CaRa + Cxa + ^ + (16) 

A simple general relativistic solution for which Cx — 
Cu = A4 = 0, Cgr = C p h, k — (flat models) reads 
as 



l/o(t) = \c vh {t-t f 



(17) 



However, it is more interesting to consider a spatially 
closed universe with k — +1 which reads as 



*(*) = 



-1 



\ C* h sm2(t-t ) 



(18) 



This immediately shows that the admissible domain for 
time is: tt/2 < (t—to) < it. It is interesting to notice that 
in (|18f) both initial (Big-Bang) and final (Big-Crunch) 
singularities emerge for infinite values of the scale factor 
(i.e., g — > 00 when a(t) — » 00). This is because g = C p ha 
in a conservation law. The new type of singularities are 
now commonly called Big-Rip singularities |20l |21| . 

It is easy to verify that the standard FRW limit can 
be obtained from 1)1(1 1) by taking 1/A — > 0, i.e., where the 
inverse of the brane tension tends to zero. The contri- 
bution from the brane p 2 term is important only when 
p > A > (lOOGeV) 4 . 

The type of cosmology we study (cf. Eq. JTJJJ) is a 
special case of a generalized brane model which is now 
commonly called a Cardassian model [29U30| |. In the Car- 
dassian model the second term on the right-hand side of 
the equation (|10l) has an arbitrary power n of the energy 
density. For brane models we study n = 2. 

In order to study observational tests we now define 
dimensionless observational density parameters |25l |26| 
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where the Hubble parameter H — a/ a, and the decel- 
eration parameter q = —da/ a 2 , so that the Friedmann 
equation l|12|) can be written down in the form 



V G R + fi A + fife + fi/ 



fi?y = 1. 



(20) 



Note that fi^ in 1)19(1 . despite standard radiation term, 
can either be positive or negative. Using 1)19(1 . the equa- 
tion ((T2^ can also be rewritten as (compare Eq.(10) of 

H) 



fiA 4 



37-2 



n G R + (3 7 - i)n A + fit 



(21) 
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It is also useful to express the curvature of spatial sections 
in terms of observational parameters by using (|2U[) and 

-n fc = ^n Gfl + 3 7 o A + 2fi w - ? - 1. (22) 

These relations <|21ll and (|22|l may allow to write down 
observational quantities for the brane models to study 
their compatibility with astronomical data which is the 
main objective of the present paper. Obviously, these 
quantities depend on the present densities of the different 
components of matter content Q given by l|19|l and their 
equations of state reflected by the value of the barotropic 
index 7. 




III. BRANE MODELS IN A 
TWO-DIMENSIONAL PHASE PLANE 

For the considerations of this section and without the 
loss of generality it is useful to simplify the equation ltTU|) 
by assuming that the parameters k 2 = 1 and 4 M 2 = 8tt 
which gives 



FIG. 1: The phase portrait of the system 1251 . There is one 
critical point - a saddle point. This critical point (xo, 0) rep- 
resents the Einstein Static Universe. The accelerating models 
lie in the domain x > xo- Therefore the acceleration region 
is situated to the right of the saddle point. The flat model 
trajectory flk,o ~ separates the regions of the models with 
negative curvature Qk,o > and positive curvature flk.o < 0. 
The line x = 1 corresponds to a = ao. 



A- 6 

«6 



so that the equation (fTTTf> reads as 



H 2 




A 5 + A 



k 2U Q 



(23) 



and it contains only three independent constant A, A5 
and Uq. The last constant comes form the conservation 
law of dark radiation (JJJ). Without loosing a generality 
in our subsequent analysis we also put A5 = 1 so that 
the equation (|23[1 takes the form 
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(24) 



As we have mentioned already brane cosmology with 
the Randall-Sundrum ansatz is mathematically equiva- 
lent to a multifluid Friedmann cosmology. It is governed 
by the following dynamical system |2^| 



y 



(25) 



dV 1 n ,„ 
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-(37*-l) 



where dot denotes differentiation with respect to a new 
rescaled time variable T:t — > T = \Ho\t; and fi^o are 
the density parameters as defined in l|19|) : x = a/do, and 
Pi = (7j — l)pj. For the sake of this Section we choose 
that a 4-dimcnsional brane universe is filled with non- 
relativistic matter (dust) on the brane only and the cos- 
mological constant, so that we effectively put 7 m = 1, 
7a = 2 (this is due to the fact that dust contribution 



which comes from non-relativistic matter scales as stiff- 
fluid in standard case - cf. eqs. i|ll |) -l(T2 )l . 7A = 0. How- 
ever, as we already mentioned, the brane (g 2 ) contribu- 
tion from other types of fluids will scale differently sim- 
ulating various types of fluids in standard cosmology. 

Here, V(x) plays the role of the potential function for 
the system (12511 



V(x) 



2 ^fc,o + ^2 ^ 4 >° 3 



-(3 7 i"2) 



(26) 



Let us note that the dynamics can also be represented 
in terms of a Hamiltonian dynamical system with the 
Hamiltonian of the form 

H = & + V(x) , 

and the system should be considered on the zero-energy 
level 



n = o. 



(27) 



The phase portrait of the system Ij25(l for the case of 
the dust 7=1 matter on the brane with flA 4 ,o = 0.69, 
^m,o = 0.3, £1a,o = 0.01 is demonstrated in Fig.Q] 

The differences in the behavior of trajectories manifest 
at high densities. Then, the structure of the dynamical 
behavior at infinity is modified. In order to illustrate 
the behavior of trajectories at infinity, the system (|25|l is 
represented in Fig. m the projective coordinates 
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(z, u) — map (28) 
(v, w) — map. (29) 
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FIG. 2: The phase portrait of the system 1251 in the pro- 
jective coordinates (z,u) for the analysis of the behavior of 
trajectories at infinity. Two critical points appear at infinity 
(0, uo±). They are of saddle type. 



FIG. 3: The dependence of Qa 4 on x. The evolutional path 
are represented by levels of constant A4. Let us note that the 
domain under the characteristic curve Qa 4 (x) is non-physical. 



The two maps (z, u): z = 0, —00 < u < 00 and (v, w):v — 
0, —00 < w < 00 cover the behavior of trajectories at the 
infinity circles x — 00 and y = 00. 
The original system 

x = P{x,y) , 
y = Q(x,y) , 

in the projective coordinates (z,u) and after the time 
reparameterization r — + t± : dr\ = xdr takes the form 

i = zP*(z, u) , 

u = Q*{z,u) — uP*{z,u) , 

where 

P*(z,u) = z 2 P(\/z,u/z) , 
Q*(z,u) = z 2 Q{\/ z,u/ z) , 

and the dot denotes differentiation with respect to time 

n- 

The representation of the dynamics as a one- 
dimensional Hamiltonian flow allows to make the clas- 
sification of possible evolution paths in the configuration 
space which is complementary to phase diagrams. It also 
makes it simpler to discuss the physical content of the 
model. Finally, the construction of the Hamiltonian may 
allow to study brane quantum cosmology in full analogy 
to what is usually done in general relativity. 

From equation (|25|l we can observe that the trajec- 
tories are integrable in quadratures. Namely, from the 
Hamiltonian constraint TL = E = we obtain the inte- 
gral 



t-tn = 



da 



a y/~2V(a) 



(30) 



Obviously, for a specific form of the potential func- 
tion H2tj|) we can obtain exact solutions. Some special 
cases have been already given in Section III. 

It is possible to make a classification of the evolution 
paths by analyzing the characteristic curves which rep- 
resent the boundary equation of the domain admissible 
for motion. For this purpose we consider the equation 
of zero velocity, d = which constitutes the boundary 
M = {a:V(a) = 0}. 

From equations (|21|l and (122[) the cosmological con- 
stant can be expressed as a function of x as follows 

^a 4 ,o(^) = -x~ 2 (n^o^ 1 + Q\,ox~ 4 + fik,o) > (31) 

where we have abbreviated Ugh by £l m . The plot of 
^A it o(x) for different k is shown in Fig. |3| Finally, we 
can consider the evolution path as a level of £1a 4 ,o = const 
and then we can classify all the models with respect to 
their quantitative properties of the dynamics. 

The present experimental estimates based on baryons 
in clusters, give £l m ,a ~ 0.3, and the location of the first 
peak in the CMB detected by Boomerang and Maxima 
suggests a flat universe. This implies that our Universe 
with the cosmological term is presently accelerating if 
(under the assumption of the dust 7 = 1 matter on the 
brane) 



The required value of f^.o seems to be unrealistic (see 
next section) and therefore the cosmological constant is 
still needed to explain the present acceleration of the Uni- 
verse if we do not accept phantom type of matter. 
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IV. A REDSHIFT-MAGNITUDE RELATION 
FOR BRANE UNIVERSES 

Cosmic distance indicators like the luminosity dis- 
tance, depend sensitively on the spatial geometry (cur- 
vature) and on the present density parameters of differ- 
ent matter components and their form of the equation of 
state. For this reason a redshift-magnitude relation for 
distant galaxies is proposed as a potential test for the 
Fricdmann brane model. 

Let us consider an observer located at r — at the 
moment t = to which receives a light ray emitted at t = t\ 
from the source of the absolute luminosity L located at 
the radial distance r±. The redshift z of the source is 
related to the scale factor a(t) at the two moments of 
evolution by 1 + z = a(to) / a(t\) = ao/a. If the apparent 
luminosity of the source as measured by the observer is I, 
then the luminosity distance 4 of the source is defined 
by the relation 

I- L L 1 (32) 

4tt4 4tt (1 + z) 2 a 2 r\ ' 1 ' 

The observed and absolute luminosities are defined in 
terms of K-corrected apparent and absolute magnitudes 
m and M. . When written in terms of m and M. , Eq. i|32|) 
yields 



m(z) =M + 5 log 10 [£>£(»], 



(33) 



where M. = M — 5 log 10 Ho + 25, and 2?l = H^d^ is the 
dimensionless luminosity distance. For homogeneous and 
isotropic Friedmann models one gets |2^| 



(34) 



where S(x) = sui\ for /C = — ttk,o', S(x) — X f° r & — 
0; S(x) = sinhx for K, — Qfc,o- From the Friedmann 
equation 112|) and the form of the FRW metric we have 



a H 



lf(z')]- 1/2 dz' 



(35) 



and the 'brany' dust contribution Q\ m .o appears in front 
of (1 + z') 6 . If there is only radiation (7 r = 4/3) on the 
brane the luminostity distance reads as 



d L( z ) = —FT- X 
tin 



dz' 



o vAp + z'Y + fi Arj0 (i + z'f + r> A4 ,o 



,(38) 



and the 'brany' radiation contribution fl\ r ,o appears in 
front of (1 + z') s . Another important example is the 
phantom matter on the brane (7 = —1/3) which gives 
the luminosity distance in the form 



d L (z) 



l + z 



H 



Jo 



dz' 



(T+I 7 ) + (T+FF + A4 '° 



(39) 



and the 'brany' phantom contribution fl\ p h,o appears in 
front of (l + z')- 2 . 

In order to compare various matter on the brane mod- 
els with the supernova data, we compute the distance 
modulus 

fi Q = 51og(d L ) + 25 

where d^ is in Megaparsecs. The goodness of a fit is 
characterized by the parameter 



E 



(40) 



fiz,i 



In il4TJ|l Hq . is the measured value, fj, i is the value cal- 
culated in the model described above, az i is the mea- 
surement error, a 2 iZ i is the dispersion in the distance 
modulus due to peculiar velocities of galaxies. 

We assume that supernovae measurements come with 
uncorrelated Gaussian errors and in this case the likeli- 
hood function L can be determined from a x 2 statistic 
£ oc exp (-x 2 /2) dill. 



where 

m 



,o(l + *') 



67 / /\ 37 

+ Qgr,q (l + z 



+fi fe ,o 1 + z) + n u . (l + z ) + ^a 4i0 



(36) 



In order to understand the variability of g 2 contribu- 
tion onto the dynamics of the brane universes which is 
formally described by the term fi^.o let us give a cou- 
ple of particular examples. For instance, taking only the 
non-relativistic matter on the brane (7 m = 1) we can see 
that the luminosity distance can be written down as 



a t \ 1+z 
d ^ z ) = —S— x 

Hq 



dz' 



lo y/n mfi (i + z') 3 + n Am , (i + Z'f + «a 4i0 



,(37) 



V. BRANE UNIVERSES TESTED BY 
SUPERNOVAE 

We first test brane models using the sample A of Perl- 
mutter SN la data. In order to avoid any possible se- 
lection effects we work with the full sample without ex- 
cluding any supernova from that sample. It means that 
our basic sample is Pcrlmuttcr sample A. We checked our 
analysis using Permutter's samples B and C but not sig- 
nificance difference was obtained. We estimate the model 
parameters using the best fit procedure. In the statisti- 
cal analysis we also use the maximum likelihood method 

m 

Firstly, we will study the case 7 = 1 (dust on the 
brane; we will label Qgr by fi m ). The case 7 = 2/3 
(cosmic strings on the brane) has recently been studied 
in [2I where, in fact, VLu and VL\ were neglected and 
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TABLE I: Results of the statistical analysis for the dust mat- 
ter (7 = 1) on the brane for the considered samples of SNIa. 
First two lines for each sample are the best fit model and the 
best fit flat model for the sample. Third line is the best fit 
flat model obtained by minimization over M, while the last 
line is the best fit flat model with fl m ,o = 0.3 obtained by 
minimization over M. 
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where the term Q m o(l + z ) 3 was introduced in order to 
admit dust matter on the brane. This case was already 
presented in a different framework in Ref . [2(| . Secondly, 
we will study the case 7 = —1/3 (phantom on the brane 
|2C| - we will label this type of matter with Q p h instead 
of Mgr)- 

Let us first estimate the value of M. from the full sam- 
ple of 60 supernovae. For the flat ACDM model we ob- 
tained M. = —3.39. We use Perlmutter data and method 
of maximum likelihood estimation on this data to esti- 



mate the different cosmological parameters of interest, 
i.e. pairs (fi m , , to\,o), (^m,o, ^A,o)- 

The result of statistical analysis is presented in some 
figures and in Table I. Fig. 0] illustrates the confidence 
level as a function of (£l m .o, ^a,o) f° r the flat model 
(fifc.o = 0) minimized over M. with Oa^o = 1 — ^m,o — 
fifc.o — r^A.o- In present cases we formally assume that 
both positive and negative values of Oa.o are mathemat- 
ically possible although the negative values are in fact 
only possible if we admit the timelike extra dimensions 
|33| . We show that the preferred intervals for f2 m) o and 
^a,o are ft m ,o < 0.4 and il\,o > 0. 

We repeat our analysis for the case of Qk,o 7^ 0. The 
results of the statistical analysis are presented in Table I. 
Two upper lines for each sample are best fit model and 
best fit flat model for sample with fixed M. as obtained 
for the flat ACDM. Third line is best fit flat model obtain 
with minimization over A4 while last line is best fit flat 
model with Q mj o = 0.3 (preffered by present extragalac- 
tic data HH) obtained with minimization over M.. 

One should note that the result obtained for Perl- 
mutter's sample C is a little bit different than the one 
obtained in preliminary analysis for M. = —3.39 
(O fe , = 0., (O m , = 0.21, fl Xfi = 0.048 ft A , = 0.75). 
It reflects the sensitivity of the results to changes of Ai 
parameter. 

The Perlmutter data were gathered four years ago, 
hence it would be interesting to use more recent super- 
novae data as well. We decided to test our model using 
this new sample of supernovae. Recently Knop et al. |35J 
have reexamined the Permutter's data with host-galaxy 
extinction correctly assessed. The mentioned authors dis- 
tinguished few subsets of supernovae from this sample. 
We consider two of them. The first is a subset of 58 su- 
pernovae with extinction correction (Knop subsample 6; 
hereafter K6) and the second one a sample of 54 super- 
novae with low extinction (Knop subsample 3; hereafter 
K3). Sample C and K3 are similarly constructed because 
both contain only low extinction supernovae. 

Another sample was presented by Tonry et al. j3|J 
who collected a large number of supernovae published 
by different authors and added eight new high redshift 
SN la. This sample of 230 Sne la was recalibrated with 
consistent zero point. Whenever it was possible, the ex- 
tinctions estimates and distance fitting were recomputed. 
However, none of the methods was able to be applied to 
all supernovae (for details see Table 8 in j3||). This sam- 
ple was improved by Barris who added 23 high redshift 
supernovae including 15 at z > 0.7 doubling the pub- 
lished number of object at this redshifts |37j . 

Despite of these problems, the analysis of our model 
using this sample seems to be interesting. We decided to 
analyze the sample of 218 SNe la (hereafter sample TBI) 
which consists of low extinction supernovae only (median 
V band extinction Ay < 0.5). 

Tonry and Barris [3(| |37j presented redshift and lu- 
minosity distance observations for their sample of super- 
novae. Therefore, Eqs. (23a) and (23b) should be modi- 
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fied H|: 

m - M = 5 log 10 (P L )To„ry - 5 log 10 65 + 25 (41) 

and 

M = -51og 10 ff + 25. (42) 

For H = 65 km s" 1 Mpc" 1 , we obtain M = 15.935. 

The results obtained for the flat model with the 
Tonry/Barris sample of 218 SNIa is very similar to that 
obtained for flat model with the Perlmutter sample. 
However, even then we allowed Q>k,o 7^ than for TBI 
sample we obtained that preffered model of the universe 
is nearly flat one, which is in agreement with CMBR data. 
It is the advantage of our model in comparision to ACDM 
model, where 0,0] the high negative value of Qk,o was 
the best fit, although £lk,a is also statistically admissible. 
In order to find the curvature, they additionally used the 
data from CMBR and extragalactic astronomy. 

We also confront our model against Knop's sample. 
While for the flat model we obtained similar results than 
for the previous samples, in the case without of any priors 
on £lk,o Knop's sample prefers its highly negative values 
together with the negative values of ^l\,o- Note, however 
that it is an unphysical case, because the branc tension 
should be positive unless timelike extra dimensions are 
admitted [33|]. One should note that with this assump- 
tion the Knop's sample suggests (il m .o < 0.3 which is in 
agreement with the result obtained for this sample in the 
case of ACDM model [35|. 

Applying the marginalization procedure over f2fc,0j 
f^A.o, M for the Perlmutter sample A we find the lowest 
value of x 2 for each pair of values (r2 m ,0j ^A,o) as shown 
in Fig. [S] The favoured intervals for (f2 m ,0)^A,o) are 
^m.o < 0.5 and f^A.o — 1-3. As the best fit we obtain: 
fi*,0 = -0.48, fl m ,o = 0, ^a,o = 0.14, and M = -3.43. 
There is a marginal difference between that result and 
the one obtained for A4 = —3.39. Another example 
is that for M = -3.39 and Q kfi = -0.20, (f2 m , = 0, 
^a.o = 0.118 Oa,o = 1-08 and for the case noted in the 
first line of tab I we obtain nearly the same \ 2 value. 

In fact, we obtained a 3D ellipsoid in a 3d parameter 
space f2 m ,07 £l\,o, ^a 4 ,o- Then, we have more freedom 
than in the case of analysis of Ref . j where there was 
only an ellipse in a 2D parameter space fl m ,o and f^A^o- 
It clearly demonstrates that the results obtained from the 
best-fit analysis should be supported by the analysis of 
the confidence levels for the parameter intervals obtained 
from maximum likelihood method. Support from CMBR 
and extragalactic astronomy results may also be useful. 

In the case of fi^o ^ 0, the confidence level for values 
of pairs (f2 m .o, ^a,o) (Perlmutter sample A) are shown in 
the standard way after minimizing over M, Qk,0i ^a,o 
in Fig. One can see that a non-zero cosmological con- 
stant is required in our model. This analisis was repeated 
with fixed value of M = —3.39 (Fig. |7J). These figures 
shows that procedure of the minimizing over M. , however 
important, is not crucial for selecting preferred intervals 
of (^ m ,Oj^A,o) m our model. 



In Fig. we presented confidence levels in the plane 
(O mj o,^A,o) for TBI sample. The allowed area is in 
agreement with that obtained for the Perlmutter sample, 
but significantly smaller. The earlier conclusion that the 
non-zero cosmological constant is required is confirmed. 
However, as we will see further, the admission of strongly 
negative pressure matter (phantom) on the brane does 
not require cosmological constant. 

Padmanabhan and Choundhury [3^, ^3 suggest that 
dividing the sample finto the low and high redshift su- 
pernovae gives interesting results. They obtained that 
although the full data sets of SNIa strongly rules out 
models without dark energy, the high and low redshift 
data sets individuall y, a dmit decelerating models with 
zero dark energy [33, • We decide to check this result 
it the case of our model. We divide TBI sample for 120 
low redshift SNIa with z < .0.25 (TB2 sample) and 98 
high redshift SNIa with z > 0.25 (sample TB3). The re- 
sults are presented in Table I and in Figs. ED- Fig.^J One 
can see that the preferred values of the model parameters 
are different in both cases. For low redshifts z < 0.25 the 
data sets taken individually results in decelerating mod- 
els without cosmological constant, while for high redshift 
z < 0.25 data sets it requires the non-zero cosmological 
constant. It confirms that pure low redshift data cannot 
be used for discrimination between cosmological models 
effectively (see [39|). 

On should note, that the present extragalactic data 
for galaxy clusters (cluster baryon fraction) with CMB 
anisotropy measurement prefer a flat model (f2fc,o) with 
f2 OT) o — 0.3 34]. For a deeper statistical analysis of brane 
expansion scenario in explaining the currently accelerat- 
ing universe we consider ID plot of the density distribu- 
tion of r^A.o- From this analysis one can obtain the limits 
at the la level. Fig.llllshows the density distribution for 
in the flat model with f2 mi o = 0.3. This distribution 
is obtained from the marginalization over M.. On the 
base of the noted above results we assumed that "true" 
value of M. is in the interval [—3.37, —3.44] and that in- 
terval we take into account during the marginalization 
procedure. One can conclude that with the probability 
of 68.3 we get fi^o < 0.019 while with the probability 
95.4 we get ft A ,o < 0.037. 

If we formally consider the possibility of a negative 
value of f2A,o which can be interpreted as the negative 
brane tension with dust on the brane p3| the correspond- 
ing distribution function is shown in Fig. 1121 One can 
see that < 0.020 on the la confidence level and 

^a.o < 0.038 on the 2a confidence level with preferred 
value f^A.o = 0.004. However, the error in estimation 
is so large that possibility that Oa,o = cannot be ex- 
cluded. One should note, that when we formally increase 
the analyzed interval of M. to [—3.29, —3.49] we obtain 
the preferred value f^.o = 0.013. Also, the errors in es- 
timation increase and J7a,o < 0.037 on the la confidence 
level and f2A,o < 0.061 on the 2a confidence level. It 
demonstrates the sensitivity of the results of estimation 
of f^A.o with respect to small changes of M. parameter. 



In the near future the SNAP mission is expected to 
observe about 2000 type la supernovae each year, over a 
period of three years ■ Therefore it could be possible 
to verify hypothesis that fi^.o > because errors in the 
estimation of SIa.o will decrease significantly. We test 
how large number of new data should influence on errors 
in estimation of CI\.q. We assume that the Universe is 
flat with Sl m>0 = 0.28, Q x ,o = 0.01 and M = -3.39. For 
the model with dust matter on the brane we generate 
the sample of 1000 supernovae randomly distributed in 
the redshift range 0.01 < z < 2. We assume a Gaussian 
distribution of uncertainties in the measured values of m 
and z. The errors in redshifts zare of order la = 0.002 
while the uncertainty in a measurement of the magnitude 
m is assumed as la = 0.15. The systematic uncertainty 



0.02 mag at z = 1.5 41] which means 
(0.02/1.5)0. For the sample generated 



limits is a sys 
that a sys (z) 

in such a way we shoul now repeat our analysis. The 
error for SIa.o — 0.0007 on the confidence level of 68.3, 
while ~ 0.0013 is on the confidence level 95.4. It is 
clearly confirmed that the error in measurement of f^A.o 
from supernovae data will decrease significantly in the 
new future. 

In Fig. El we present the plot of residuals of redshift- 
magnitude relationship for the supernovae data for 
^m,o = 0.3. We present the residuals between the 
Einstein-de Sitter model and Perlmutter model, best fit- 
ted flat brane model and best fitted brane model. As a 
result, with the increasing impact of f^A.o (higher Oa,o) 
the high-redshift supernovae should be brighter than the 
expected by the Perlmutter model. Let us note that for 
the best-fit value of 11a, = 0.020, the difference between 
the brane model and the Perlmutter model should be de- 
tectable for z > 1.2 (because Am > 0.2). One should 
observe that there is a small difference between the best 
fit flat model (with minimal value of x 2 ) where we obtain 
^a.o = 0.006 and the value obtained from the minimiza- 
tion procedure fiA.o = 0.004 (see Fig. lT^ . 

For comparison we presented the best fits (Fig. IT5|> 
without a specifically assumed value of £l m ,o- This lat- 
ter case is, characterized by greater expected differences 
between the brane model and the Perlmutter model for 
high redshift supernovae. 

The above prediction for brane models that the high 
redshift supernovae are brighter than expected in ACDM 
(Perlmutter) model could also be detectable by the new 
SNAP data. It gives a possibility to discriminate between 
predictions of the two models. Of course similar effect is 
expected for galaxies, but for an extended object such an 
effect is much more difficult to detect than for the point 
sources. 

In Fig. El we present the difference in redshift- 
magnitude relation ll-iol) for brane models with phantom 
matter on the brane (7 = —1/3) using Perlmutter's sam- 
ple A. Note that the theoretical curves are very close 
to that of 0| which means that the phantom cancels 
the positive-pressure influence of the g 2 term and can 
mimic the negative-pressure influence of the cosmological 
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FIG. 4: Confidence levels in the plane (f2 m ,o, Oa.o) minimized 
over M for the flat model, and with Qa 4 ,o — 1 — fim.o — fifc.o — 
f2 A ,o- The figure shows the ellipses of the preferred values of 
fi m ,o and Qa 4 ,o- The results prefer positive values of fiA,o, 
but negative values are allowed as well (Perlmutter sample 
A). 
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FIG. 5: The levels of constant x 2 m the plane (f2 m ,o, Qa 4 ,o) 
marginalized over A4, and with f2A 4 ,o = 1 — O m ,o — ^*,o — 
Q,\ t o The figure shows the preferred values of Qm,o and Qa 4 ,o 
(Perlmutter sample A). 



constant to cause cosmic acceleration. From the formal 
point of view the best fit is (x 2 = 95.4) for £lk,o = 0.2, 

n ph , = o.7, n A ,o = -o.i, n u = 0.2, n Aifi = which 

means that the cosmological constant must necessarily 
vanish. From this result we can conclude that phantom 
matter p = — (4/3)p can mimic the contribution from 
the A 4 -term in standard models. For the best-fit flat 
model (flk.o — 0) we have (\ 2 — 95.4): Q p h.o = 0.2, 

ftA, = -o.i, n u = 0.2, r*A 4 , = 0.7. 
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FIG. 6: Confidence levels in the plane (f2 m ,o, ^a 4 ,o) mini- 
mized over Ai, and with fiA 4 ,o = 1 — fim.o — flk.o — ^a,o- The 
figure shows the ellipses of the preferred values of f2 m ,o and 
^a 4 ,o (Perlmutter sample A). 



FIG. 7: Confidence levels in the plane (flm.o, flA 4 ,o), M = 
—3.39, and with Qa 4 ,o = 1 — flm,o — fifc.o — Oa,o. The figure 
shows the ellipses of the preferred values of fi m ,o and Oa 4 ,o 
(Perlmutter sample A). 



VI. ANGULAR DIAMETER SIZE AND THE 
AGE OF THE UNIVERSE 

Another cosmological test which can be performed is 
the angular diameter size 

9=«zpr (43) 

of a galaxy with a linear size d. In a flat dust (7 = 1) 
universe 9 has the minimum value z m i n = 5/4. It is par- 
ticularly interesting to notice that for flat brane models 
with fl\ « 0, flA 4 ~ the dark radiation can enlarge the 
minimum value of 9 while the ordinary radiation lowers 
this value (25|, i.e., 

Zmin = ^{ n u-l + V'^u + l) > ~ (44) 

for flu < 0. This is a general influence of negative dark 
radiation onto the angular diameter size for brane mod- 
els. One can also notice that there exists a restriction on 
the amount of negative dark radiation coming from i|44|) 
(^w > — 1/3) which can serve as a test for the admissible 
value of flu = —1/3 (z m i„ = 2) in order to observe the 
minimum. 



More detailed analytic and numerical studies show that 
the increase of z m i n is even more sensitive to negative 
values of fl\ which, unfortunately, are admissible only for 
timelike extra dimensions |3^| . Similarly as for the dark 
radiation flu, the minimum disappears for some large 
negative fl\. Positive flu and fl\ make z m i n decrease. 
In Fig. El we present a plot from which one can see the 
sensitivity of z m i n to flu- 

As far as a possible contribution of fl\ is concerened, 
then from Fig. ^| we can see that its influence on the 
angular size 9 is relatively weak. As a conclusion we can 
say that it is possible to test values of flu,o and fl\o from 
the angular diameter minimum value test, but because of 
the evolutionary effects and observational difficulties the 
predicted differences are too small to be detect. We have 
also checked that phantom matter fl p h has very little 
influence onto the value of z m i n . 

Now let us briefly discuss the effect of brane parameters 
and phantom matter onto the age of the universe which 
according to ifT^I is given by 

#0*0= / {fl G RflX^ +i + fl X , Q X-^ +i 

Jo 

+fl u .o + MkfiX 2 + r2A 4 ,o£ 4 }~" xdx, (45) 
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FIG. 8: Confidence levels in the plane (Q m ,o, ^a 4 ,o) for the 
sample Tonry/Barris 218 low extincted SNIA. 



where x — a/a . We made a plot for the dust 7 = 1 
on the brane in Fig^| which shows that the effect of 
quadratic in energy density term represented by tt\ is to 
lower significantly the age of the universe. 

For the dust matter and the cosmological constant the 
age of universe is given by 



to = — — In v (46) 

3-ffo\/^A,0 "m,0 + \/4fi.\,0^A,0 

In Fig.|2]we made another plot of (|46|l in Gyrs for the fiat 
model for different We can again see that the dust 

matter on the brane which mimics brane contribution as 
stiff-fluid decreases significantly the age of the universe. 

The problem can be avoided, if we accept phantom 
7 = —1/3 on the brane |2jj, since the phantom has a 
very strong influence to increase the age. In Fig. EU| 
we made a plot for this case which shows how phantom 
energy enlarges the age. 

Therefore the cosmological constant or the phantom 
matter is needed to explain the problem of the age of 
the universe. 



FIG. 9: Confidence levels in the plane (Q m ,o, ^a 4 ,o) for the 
sample Tonry/Barris 120 low extincted SNIA with z < 0.25. 
The figure illustrate that with sample of low redshift super- 
novae we obtain that model without cosmological constant is 
statistically admissible. 

VII. DOPPLER PEAKS 

The CMB peaks arise from acoustic oscillations of the 
primeval plasma. Physically, these oscillations are rep- 
resented by hot and cold spots in CMB temperature. A 
wave which had a maximum of the density at the time of 
the last scattering corresponds to a peak in the power 
spectrum. In the Legendre multipole space this cor- 
responds to an appropriate angle 9 a subtended by the 
sound horizon at the last scattering surface. Higher har- 
monics of the principal oscillations, which have oscillated 
more than once, correspond to secondary peaks. 

For our purpose it is very important that the location 
of these peaks are very sensitive to the variations in the 
model parameters. Therefore, it can be used as a probe to 
constrain the cosmological parameters and discriminate 
among various models. 

If 9 a is the angular size of the density fluctuation, then 
the acoustic scale for the peaks is [l a = ir/9 a ) [3l| 

Z f[f(z)]-^dz 

la = n^ (47) 

/ c s [f{z'))-VHz' 

Zdec 
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FIG. 10: Confidence levels in the plane (f2 m ,o, ^A 4 ,o) for the 
sample Tonry/Barris of 98 low extincted SNIa with z > 0.25. 

{zdec — 1100) where c s is the speed of sound and the 
position of the peaks is 



k = l a {i ~ 4>i) 
(i = 1,2,3), where the phase shift [HE 



0.267 



r(zd cc ) 



-i o.i 



0.3 



(48) 
(49) 



The phase shift is caused by the pre-recombination 
physics (plasma driving effect) and, hence, has no sig- 
nificant contribution from the term containing brane in 
that epoch so that it was taken from standard cosmology. 
We take fi b ,o^ 2 = 0.02, r(z doc ) = Pr(z dcc )/ p m (z dcc ) = 

^r,o(l + z dcc)/^m,0, ^r,0 = ^7,0 + ^7,0 = 2.48/l 2 • 

10~ 5 , flvfi = 1.7h~ 2 ■ 10 -5 , r(zdcc) is the ratio of radia- 
tion to matter densities at the surface of last scattering 
and we have O ri o = 9.89 • 10~ 5 , Slb.o = 0.05, and the 
spectral index for initial density perturbations n = 1, 
and h = 0.65 Q. 

For dust matter f2 m .o, radiation £l r ,o, and dark radia- 
tion £lu,o (^7,0 = ^r.o + Qu,o) on the brane under the 
assumption that we neglect the influence of the quadratic 
in energy density term of radiation f2 rj o we have the ef- 
fective speed of sound 



-eff 



= c 2 + 



n 6 ,o(i + z) 3 + |n 7 , (i + z) 4 
n Xo + n 6 , (i + zf + n 7 , (i + zf ' 



(50) 
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FIG. 11: The density distribution for in the model with 
(only positive) brane fluid present. We obtain the limit f2x,o < 
0.019 at the confidence level 68.3, and fiA,o < 0.037 at the 
confidence level 95.4 (fi m ,o = 0.3, Perlmutter sample A) 




FIG. 12: The density distribution for Q\ t o in the brane model. 
We obtain that fi A , < 0.020 at la level and fi A , < 0.038 
at 2a level. Both positive and negative values of Q\,o are 
formally possible. For simplicity, we only mark limits on the 
positive side. (fi m ,o = 0.3, Perlmutter sample A). 



with c 2 of the same form as in standard cosmology and 



O 2 

A '° 2f2 A ,o 



(51) 



The influence of brane on the location of the peaks is 
to shift them towards higher values of I. For example, for 
^m,o = 0.3, Ob,o = 0.05, h = 0.65, the different choices 
of il\.o yield the following 

^A,0 = 0: t pea k,l = 225, Zpoak,2 = 535, Zpeak,3 = 847 

^A,0 = 1-5 • 10 lj : ^poak.l = 227, Z pca k,2 = 540, ? pe ak,3 = 853, 
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FIG. 13: The density distribution for Qx,o in the brane model 
(with M in the interval [-3.29; 3.49]). We obtain that fi A ,o = 
0.013 and fi A ,o < 0.037 at la level and n x ,o < 0.061 at 2cr 
level. Both positive and negative values of S1a,o are formally 
possible. For simplicity, we only mark limits on the positive 
side. (fini,o = 0.3, Perlmutter sample A). 
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FIG. 14: Residuals (in mag) between the Einstein-de Sitter 
model and: the Einstein-de Sitter model itself (zero line), 
the Perlmutter flat model (upper curve), the f2 m ,o = 0.3 the 
best-fit flat brane model (fi A ,o = 0.006 and Qk,o = 0) (upper- 
middle curve), and fi m ,o = 0.3 best-fit brane model (fi^.o — 
0.020, flk,o = —0.1 (lower- middle curve). 



We could also analyze the influence of dark radiation 
term in the brane world cosmology. The corresponding 
term in this case scales just like radiation with a constant 
Po and both positive and negative values of p r ,o (pu,o) 
are possible. Dark radiation should strongly affect both 
the BBN and CMB. Ichiki et al. used such obser- 
vations to constrain both the magnitude and the sign 
of dark radiation in the case when the p 2 term coming 
from the brane is negligible. We take their limits on 
dark radiation with negative contribution because of the 
tension between the observed 4 He and D abundance [ifij . 
In our case we take into consideration the brane term p 2 
and obtain the following positions of the first three peaks: 
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FIG. 15: Residuals (in mag) between the Einstein-de Sitter 
model and: the Einstein-de Sitter itself (zero line), the Perl- 
mutter flat model (upper curve), the best-fit flat brane model 
(upper-middle curve), flk,o = 0, f2 m ,o = 0.25, Cl\ t o = 0.02, 
^a 4 ,o = 0.73 and the best-fit brane (lower- middle curve) 
fifc.o = -0.9, On.o = 0.59, fi A ,o = 0.04 ft A4 ,o = 1-27. 
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FIG. 16: Residuals (in mag) in redshift-magnitude relation 
for 7 = —1/3 brane universes (phantom matter on the brane). 
The top line is the best-fit flat model of Ref. [lfll ] with Q. m ,o = 
0.28, f2A (4 ) ,o = 0.72. The bottom line is a pure flat model with 
Oa (4) ,o = 0. Between - two brane models with fl\,o ^ 0: lower 
- the best-fit model; higher - the best-fit flat model. The brane 
phantom plots are very close to the top line of Ref. 0] . 

^u,o = — 1.23Q 7) o f2x,o = 10~ 12 : / p0 ak,i = 208, ? pca k,2 = 495, 

/peak, 3 = 781, 

^W,0 = —^7,0 ^A,0 = 10~ 12 : i P eak,l = 214, i pea k,2 = 545, 
/peak, 3 = 861, 

ttu.o = -0.41fi 7>0 r>A,o = 10- 12 : /pcak,i = 229 Z pcaM = 545, 

/pcak,3 = 861. 

From Boomerang observations 0] we obtain Z poa k,i = 
200 — 223, /peak. 2 = 509 — 561. Therefore, uncertainties in 
values /peak can be used in constraining cosmology with 
brane effect, namely 

^a,o < 1.0- 10~ 12 

from the location of the first peak. 

We also compare the results from the above proce- 
dure with recent bounds on the location of the first two 
peaks obtained by WMAP collaboration [47l l4e| namely 
Z peak)1 = 220.1 ± 0.8, Z P cak,2 = 546 ± 10, together with 
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FIG. 17: The angular diameter 9 for ft A = 0.1,fi ro = 
0.3, Qa 4 = 0.72, and the two values of Q.u = 0.1, —0.1 (top, 
middle) in comparison with the model of Ref. fig| with 
f2 m = 0.28, S1a 4 = 0.72 (bottom). 
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FIG. 18: The angular diameter 9 for the flat brane model for 
^m.o = 0.3 and £l\,o = 0.1,0.02,0 (top, middle, bottom, re- 
spectively). The minima for these cases are 1.04, 1.30, 1.605, 
respectively. The brane fluid causes the minimum to move left 
(towards to lower z) and the minimum value of 9 increases. 

the bound on the location of the third peak obtained by 
Boomerang collaborations Z poa k,3 = 825^^ which lead to 
quite strong constraints on the model parameters. These 
constraints can be summarized as follows. If we assume 
no dark radiation, the brane model is in agreement with 
observations provided we take 

^a,o < l-O- 1(T 15 

We have made similar considerations for phantom on 
the brane (with matter and radiation included - oth- 
erwise the speed of sound becomes imaginary, compare 
Ref. [2g|) and obtained the limit fi Aj0 < 1.0 x 10" 13 for 
Cl m = 0.3, Qph,o = 0.7 — f2 7 ,o, ^fc.o — and the location 
of the peaks is at l\ = 225, 12 = 535, ^3 = 845. 

However, the phase shift <j) in the above considerations 
is taken from the standard cosmology, i.e., we assume 
that the contribution from the brane term is insignificant 
at the pre-recombination epoch. If this assumption was 
not valid then the limit from CMB would change. 

VIII. BIG-BANG NUCLEOSYNTHESIS 

It is clear from the Friedmann equation l|10|l that the 
brane models with positive pressure 7 > 1 matter lead to 



2.50 




0.00 1 

0.00 0.20 0.40 0.60 0.S0 1.00 

Q 

FIG. 19: The age of the universe to in units of H^ 1 for the 
brane models with dust (0 < f2 m ,o < 1 on the horizontal 
axis). Here Qu,o = ^fc.o = 0, Q\,o = 0,0.05,0.1 (top, middle, 
bottom). 
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FIG. 20: The age of the universe to in units of Hq 1 for the 
brane models with phantom on the brane (0 < £l p h.o < 1 on 
the horizontal axis). Here Qu,o = 0.2, Qa,o = 0.05,0 (top, 
bottom) which shows weaker influence of the brane effects to 
increase the age. 



a dominance of the g 2 term in the early universe. How- 
ever, its admittance, even for the late universe can be 
dominant. In particular, if the restriction we found for 
dust matter on the brane from supernovae ~ 0.01 
is to be applied, then the brane term can dominate the 
universe already at redshift z ~ 2. 

For large z the brane term which comes from radiation 
scales like (1 + z) 8 and it dominates over the standard 
radiation term. This brings a potential trouble since in 
such a model radiation domination never occurs and all 
BBN predictions fail. In other words, the preferred value 
obtained from SN la data gives the p 2 term which is 
far too large to be compatible with BBN, if we assume 
that brane models do not change the physics in the pre- 
recombination epoch. 

The consistency with BBN seems to be a crucial is- 
sue in brane cosmology 0,153, El- For this reason, we 
should admit that the contribution of brane fluid fl\.o 
cannot dominate over the standard radiation and dark 
radiation terms before the onset of BBN, i.e., for z = 10 8 
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FIG. 21: The age of the universe to in units of 10 9 yrs for 
the flat brane model with fi A ,o = 0, 0.004, 0.020 (top, middle, 
bottom). 

which gives a condition 
ft 2 

n >,o 05^(1 + z ) 8 < («r,o + nw,o)(l + (52) 

where fi Wj0 < 0.11fi r)0 E3- 
Therefore, the term 

n Xfi ^(l + zf, (53) 

which describes brane effects is constrained by the BBN 
because it requires a change of the expansion rate due 
to this term to be sufficiently small, so that an accept- 
able hclium-4 abundance is produced. Taking this into 
account we obtain the following limit in this case 

f^A,o < 1-0 • 10~ 27 if z~10 8 . 

However, the situation can significantly be changed if 
we admit phantom matter in the universe because the 
phantom brane term never dominates at large redshifts. 
Quite the contrary, it eventually dominates present and 
future evolution and causes no change of BBN predic- 
tions. This is because it effectively leads to a reversal of 
eq. (152(1 in a similar way as it reverses the luminosity 
distance relation in eq. 1)39(1. 

IX. CONCLUSIONS 

We have investigated observational constraints for the 
Randall-Sundrum scenario of brane world cosmology. 
The motivation to study new observational constraints 
for some particular scenarios of the brane world cosmol- 
ogy is the demonstration that such models inspired by 
recent developments in the particle physics can be tested 
by astronomical observations. 

We have shown that this scenario is compatible with 
the most recent observations of SN la. Moreover, we 
demonstrate that only a new (high z) and more precise 
set of observations can show whether the considered class 



of models constitute a viable possibility for the descrip- 
tion of the present acceleration of the Universe. 

The brane model with dust 7 = 1 matter on the 
brane fits the supernovae la data very well (including 
SN 1997ff at z ~ 1.7). However, the cosmological con- 
stant is still required. We obtain the best-fit non-flat 
model with 17a j0 — 0.01, fi^o — —0.9, £l m .o = 0.6, 
51a, = 1-3. For the flat model with f2 m .o = 0.3 we obtain 
n A ' = 0.004 ± 0.016 as a best-fit. Whereas the best-fit 
non-fiat model is not realistic (because of large negative 
curvature), the flat model with the brane for the realistic 
value of fi m = 0.3 is in agreement with the SN la data. 

On the other hand, brane models with phantom 7 = 
— 1/3 matter on the brane also fit supernovae data and 
they may mimic the contribution from the cosmological 
constant which is then not required. 

We have demonstrated how the values of the estimated 
parameters depend on the division of the sample on high 
and low redshifts subsamplcs. As a results we obtained 
that although the full data set (Tonry/Barris sample) 
of SNIa require the cosmological term, the low redshifts 
z < 0.25 data set, admits decelerating model without 
cosmological constant. We have also demonstrated the 
sensitivity of results with respect to small changes of A4 
parameter. 

It is interesting that brane models with dust 7 = 1 
matter for high redshifts predict brighter galaxies than 
the Perlmutter model. Therefore, the difference between 
the Perlmutter model and the brane model may be de- 
tectable for high redshit z > 1.2 supernovae. 

Let us note that present data suggest that Q^qLO ~ 
10~ 2 but because of the large error of this estimation it 
is possible that f^A.o — 0. However, from future observa- 
tions of supernovae at high redshift we can expect that 
errors in estimation of Q\,o become significantly smaller. 
At present the brane theory can be neither confirmed 
nor ruled out. In this way the future SN la data should 
allow to verify the hypothesis that is so large as 
^a.o-O.OI. 

We also found the other limits on the value of fi>,,o 
from the measurements of CMB anisotropies and BBN. 
We obtain the strongest limits in these cases, namely 
fl x , < 1.0 • 10~ 12 from CMB and Q x ,o < 1.0 ■ 10~ 27 
from BBN. However, let us note that because the errors 
in estimation of Oa.o from supernovae are so large the 
results obtained from SN la do not contradict those of 
BBN and CMB. 

Of course BBN as well as CMB are very well tested ar- 
eas of cosmology which do not allow for a significant devi- 
ation from the standard model and standard expansion 
law, except at very early times. Although consistency 
with BBN and CMB is a crucial issue in brane models, 
we must remember that in such an approach we assume 
that brane models does not change the physics in the pre- 
recombination epochs. On the other hand, the results 
of analysis SNIa data are independent of the physical 
processes in the early universe. Therefore weaker lim- 
its obtained from SNIa observations may even be more 
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valuable. 
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